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ABSTRACT 


This Progress Report covering the period of June 1, 1993 to December 
1 , 1993 presents the development of an analytical solution to the heavy 
ion transport eguation in terms of a one— layer Green s function formalism. 
The mathematical developments are recasted into an efficient computer 
code for space applications. The efficency of this algorithm is accom- 
plished by a nonperturbative technique of extending the Green's function 
over the solution domain. The code may also be applied to accelerator 
boundary conditions to allow code validation in laboratoty experiments. 
Results from the isotopic version of the code with BO isotopes present 
for a single layer target material, for the case of an Iron beam 
projectile at 600 MeV/nucleon in water is presented. 


INTRODUCTION 


Future NASA missions will be limited by exposure to space radiations 
unless adequate shielding is provided to protect men and equipments 
from such radiations. Adequate methods required to estimate the damage 
caused by such radiations behind various shields can be evaluated prior 
to commitment to such missions. 

From the inception of the Langley Research Center heavy ion (HZE) 
shielding program (refs. 1-3), there has been a continued, close relation- 
ship between code development and laboratory experiment (ref. 3). Indeed, 
the current research goal is to provide computationally efficient high 
charge and energy ion (HZE) transport codes which can be validated with 
laboratory experiments and subsquently applied to space engineering design. 
In practice, two streams of code development have prevailed due to the 
strong energy dependence of necessary atomic/molecular cross sections and 
the near singular nature of the laboratory beam boundary conditions (refs. 
4-6). The atomic/molecular cross section dependence is adequately dealt 
with by using the methods of Wilson and Lamkin (ref. 7), allowing effici- 
ent numerical procedures to be developed for space radiations (refs. 

6 , B-10 ) . Although these codes could conceivably be applied to the labo- 
ratory validation, methods to control truncation and discretization errors 
would bear little resemblance to the space radiation codes attempting to be 
validated. Clearly, a radical reorientation is required to achieve the 
validation goals of the current NASA space radiation shielding program, 
and such an approach is the main thrust of this research and is briefly 


described below. 



A useful technique in space radiation shielding is the use of the impulse 
response Green's function (refs. 11,12) , which satisfies the Boltzman 
equation of the form 
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where Gj m reduces to a monoenergetic unidirectional function at the bound- 
ary, Sj(E) is the stopping power, is the total cross section, and tfjk 
is the inclusive differential cross section. An arbitrary solution to the 
Boltzman equation within a closed convex region can be written as 
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where f ffl (E',n',r) is the incident flux at the boundary (ref. 11). Since 
transport problem is formulated in terms of a single Green s function 
algorithm, the validation of the Green's function in the laboratory 
meets the objective of having a space validated code. Since there is 
hope of a Green's function based on an analytical solution of the 
Boltzmann equation (ref. 13), the resulting evaluation of the shield 
properties should be computationally efficient. 


The first step in this process is to develop an equivalent Green's 
function algorithm in one dimension to match the current capability 
in space radiation transport calculation (refs. 6,14). The algorithm 
is based on the closed form solution to the one dimensional equation 
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for a monoenergetic beam at the boundary. The probality of validation 
for 2 °Ne beams of this algorithm (with multiple scattering corrections) 
has already shown good correlation (refs. 5,15), but improvements in 
the nuclear data base are required for achieving higher correlations 
with experiment. If considerations are restricted to multiple charged 
ions then the right hand side of equation (3) can be further reduced to 
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for which a solution is presented below. 


APPROXIMATE GREENS ' S FUNCTION 

Equation (4) can be simplified by transforming the energy into the 
the residual range as 
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and defining new field variables as 
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bo that equation (4) becomes 



a a 

dx drj ^ 


^'m( I i r ji r m) 


|/, a jk&km ( z » r j i r m ) 

Ifc * 

with boundary condition 

^jm(®i r j» r m) = ^jm ^( r j ~ r m) 
and 


( 8 ) 


( 9 ) 


ipj(x,rj) = / ^im( :e * r i» r m)/m(rm)J»'!n 0°) 

m •' 0 

The solution to equation (8) may be written as 
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where zeroth order term of equation (11) is 
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and the first order term of equation (11) is 
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The second order terms of equation (11) are 
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with the condition that C^|(x, r jt r' m ) are nonzero for 
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The third order terms of equation (11) are 
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and similarly for higher order terms. In the above the g's of n arguments 
are given by 
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and 
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In terms of above , the s 


olution to equation (4) may be written as 
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for i=l, the denominator of equation (23) is 
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and for i> 1 , the denominator becomes 
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in equation (22), F„(E) i. the integral flux 
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Implementation of equation (22, can nee he accomplished independent 
of the character of the boundary values f m (<) and will give accurate 
results for both space and laboratory applications. 


DISCUSSION OF RESULTS 


Values of collision related fluxes for 3 depths of water target 
for a mono-energetic beam of Iron projectile with 600 MeV/nucleon 
are shown in figures 1 through 3 for both the perturbative and 
nonperturbative Green's function methods for the direct comparison 
of the two methods. The first three collision terms and the sum of 
all collision terms of both theories at a depth of 5 cm of water 
are shown in figure l.A through l.H. The differences in the spectral 
shape is due to the simplification of the attenuation term in the 
nonperturbative theory. The nonperturbation terms represnt the 
average spectrum while perturbation theory retains the spectral shape. 
Figures 2. A through 2.H and 3. A through 3.H are the corresponding 
comparison of the two methods at depths of 10 and 15 cm of water. 

Direct comparsions of figures 1 through 3 shows that the sequence Of 
perturbation terms appear to be converging to a result similar to that of 
nonperturbative result . 

The main advantage of nonperturbative methods are m their 
computational efficencies. The computational time required for the 
nonperturbative code is about 10 minutes on VAX 4000 compared 
to 15, 45, 90 minutes for the l-st, 2 nd and 3-rd collision terms 

of the perturbation solution. 

Figures 4. A through 4.C show the corresponding differential LET 
spectrum using the method of reference 16. The highest LET peak is due 
to the primary beam and the ion fragments. The successive peaks below 
iron are due to lower atomic weight fragments. Such LET spectra can 
be compared to experimental measurments directly. 
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Figure l.A. 1-st term nonpertubation solution at 
for a 600 MeV/nucleon Iron projectil 
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Figure l.E. 3-rd term nonperturbation solution at a depth of 5 cm of water 
for a 600 MeV/nucleon Iron projectile. 





Figure l.F. 3-rd term perturbation solution at a depth of 5 cm of water 
for a 600 MeV/nucleon Iron projectile. 
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Figure l.G. All terms nonperturbation solution at a depth of 5 cm of water 
for a 600 MeV/nucleon Iron projectile. 
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Figure l.H. All terms perturbation solution at a depth of 5 cm of water 
for a 600 MeV/nucleon Iron projectile. 
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1 Figure 3. A. 1-st term nonperturbation solution for a depth of 15 cm of water 

, for a600 MeV/nucleon Iron projectile. 
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Figure 3.F. 3-rd term perturbation solution for a depth of 15 cm of water 
for a 600 MeV/nucleon Iron projectile. 
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Figure 3.G. All terms nonperturbation solution for a depth of 15 cm of water 
for a 600 MeV/nucleon Iron projectile. 
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Figure 4.B. Differential LET spectrum for a depth of 10 cm of water for a 600 
MeV/nucleon Iron projectile. 
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Figure 4.C. Differential LET spectrum for a depth of 15 cm 
MeV/nucleon Iron projectile. 
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